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Abstract 

We describe gravitational lensing by a gravitational wave, in the regime in 
which multiple images of a light source are created. We adapt the vector for- 
malism employed for ordinary gravitational lenses to the case of a non-stationary 
spacetime, and we derive an approximate condition for multiple imaging. It is 
shown that certain astrophysical sources of gravitational waves satisfy this condi- 
tion. 
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We study gravitational lensing with a gravitational wave acting as the lens. It is 
well-known that exact gravitational waves deflect light, and light propagation through 
linearized (cosmological) gravitational waves outside the laboratory has also been studied 
IJ, expecially in conjunction with their frequency shift effect on the photons of the cosmic 
microwave background. It has also been suggested that gravitational waves may create 
multiple images of distant light sources @. We concentrate on the latter aspect of lensing 
by gravitational waves, for two reasons: i) the possibility of multiple images is associated 
to strong amplifications of the light source, which makes easier to detect gravitational 
waves; ii) a previous analysis |§ concluded that the amplification is negligible. This 
conclusion was based on the Raychaudhuri's equation, and it is not valid when multiple 
images are involved. 



1 The vector formalism for lensing gravitational waves 

We consider a gravitational wave localized in a region of space between a light source 
and an observer. The spacetime metric is g^ v = r/^ + h^ v , where rj^ is the Minkowski 
metric and {h^l <C 1. Let us consider a light ray whose unperturbed path is parallel 
to the z-aads, with 4-momentum p^ = p^ + 5p^ = (1 + 5p° , Sp 1 , 5p 2 , 1 + Sp 3 ), where 
Sp^ are small deflections. We work in the geometric optics approximation, which holds 
if \.w. >> K.m. and X e . m , > X g , w 2 /D L (where X gM . and A em . are the wavelengths 
of the gravitational wave and of the photon, respectively, and Dl is the observer- lens 
distance). We set the geometry as customary in gravitational lens theory, using lens and 
source planes with coordinates (x, y) and (s x , s y ) respectively. The action of the lens 
is described by the plane-to-plane mapping x A i — ► s A (A = 1,2), with s given by the 
"lens equation" 

s A = x A + 6p A (x) (1.1) 

(Dls and D$ are the lens-source and the observer-source distances respectively). 
The map described by Eq. ( |1.1[) has the Jacobian matrix 

Js\_ (ds^\_(l + Dd x (5p*) Dd y (5p*) \ 
\x) \dx B ) \ Dd x (5pv) 1 + Dd y (6py))> [L - A) 

where D = DlDls/ Ds- The inverse matrix A = J -1 represents the amplification 
tensor, while its determinant A = Det(J) -1 is the (scalar) amplification. Since the 
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surface brightness is conserved during lensing whenever geometric optics holds and 
^ area of an infinitesimal region in the lens plane 

area of the corresponding region in the source plane ' 

A has also the meaning of the ratio of light intensities with and without the lens. A 
small circular source will be imaged into a small ellipse whose eccentricity e is given by 
the ratio of the eigenvalues e± of A: 



e 2U/2 



(1.4) 



The vanishing of the Jacobian Det( J) indicates the failure of invertibility of the map ([TT 
Therefore, the occurrence of multiple images is signalled by the vanishing ofDet(J), and 
this is the condition that we study in the following. 

The Jacobian determinant Det(J) can be written as follows: 

Det(J) = 1 + 4J{a) D s + f{a)D 2 s \d x {5p x ) • - d y {8p x ) ■ d x (6p v )} = 

= 1 + Ji + J 2 , (1.5) 

where a = Dls/D$ and f(a) = a 2 (l — a) 2 is symmetric about a = 1/2, where it 
assumes its maximum value 1/16. The deflection Sp 1 *, computed from the equation of 
null geodesies, is 

1 /-Observer 

5p A = - / dz (h 00 + 2/ios + M' A + 0(h 2 ) . (1.6) 

2 J Source 

Moreover, 8a(Sp a ) = to first order and J\ ~ h 2 D/P < J 2 ~ h 2 {D/P) 2 for large 
values of D/P. Thus, in order to have multiple imaging, it must be J 2 < and 

r n 2 

f{ct) Ds 0a{Sp b ) ~ 1. For ordinary gravitational lenses, the probability of lensing 
of a distant source is maximum when the lens is halfway between the source and the 
observer, hence we take f(a) in the range - ^, near its maximum. Then, in order 
to have multiple imaging, it must be 

p > Sc , (1-7) 

where P is the period of the gravitational wave, S c = (4 — 10) c/D. The approxi- 
mate condition for multiple imaging ( |1.7| ) involves the "strength" h and the "size" P of 
gravitational waves, and the geometry of the problem (through D). ( |1.7|) is somewhat 
analogous to the well-known condition for multiple imaging by ordinary gravitational 
lenses, E > S c = ^ ■ j^- . 
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2 Lensing by gravitational waves and the Fermat 
principle 

The vector formalism for ordinary gravitational lenses requires the stationarity of the 
lens potential; clearly this hypothesis is not satisfied in our case, and a new approach is 
needed. We substitute the true photon path in the 3-dimensional space with a zig-zag 
path composed of two straight lines from the source to the lens, and from the lens to the 
observer. By applying a new version of the Fermat principle valid for non(conformally)- 
stationary spacetimes Q, it was shown in Ref. || that the zig-zag paths are extrema of 
the travel time functional 



1 

t = constant H 

2c 



D L D 



L^LS 



Dq 2 f° 

{x-s) + dz (h 00 + 2h 03 + h 33) 
Js 



(2.1) 



The lens equation (|1 . 1| ) and the deflection ( |1.6|) are obtained by requiring stationarity of 



this functional: V^t = 0. Therefore, the conclusions of the previous section are justified 
in the context of a rigorous formalism. 



3 Comparison between a gravitational wave and an 
ordinary gravitational lens 

We compare the action of a gravitational wave with that of an ordinary gravitational 
lens. The latter is a mass distribution described by a Newtonian potential $ (satisfying 
the Poisson equation V 2 $ = 47rp, where p is the lens mass density). The plane-to-plane 
map describing the lens action is given by the lens equation and the Jacobian matrix 
can be written as 

J=f 1-X .: A , -f.A, (3-D 



with 



X = (3 - 2) 

s *L*w-w) ' (3 - 3) 

D r+°° 9 2 $ 
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where D and S c have been defined in the previous sections and 



dip. (3.5) 

J — oo 

The Jacobian determinant is given by 

Det(J) = (l- X ) 2 -(A 2 + /x 2 ). (3.6) 

The convergence x describes the action of matter, while A and \i describe the action of 
shear. For a lensing gravitational wave one obtains 

where, to first order, one has 



Ai = 


1) 

~ ~2 


r° 

/ d\ h xx ' x , 
Js 


(3.8) 


A 2 EE 


D 
~ ~2 


f° 

/ dX h\\ t y , 
J s 


(3.9) 


Mi = 


D 
~ ~2 


dX h\\' x >y , 


(3.10) 


^2 


D 

" 2~. 


r° 

J s dX h xx ' y >x , 


(3.11) 



where the integrals are computed along the photon's path from the source to the observer. 
To first order, A 2 = — A x = — A gw and /ii = /i 2 = p gw . To the lowest order, we have 



j - = ( i +X. ) • (3 - 12) 



The convergence term is absent and hence the lens action is due only to the shear. This 
result was derived in Ref. using the Raychaudhuri equation and the optical scalars 
formalism, and is now recovered in the vector formalism. 

The deflection angle does not depend on the frequency of the light: gravitational 
waves are achromatic lenses, like ordinary gravitational lenses. However, while the lat- 
ters do not shift the frequency of the photons propagating through them, lensing gravi- 
tational waves do. In addition, gravitational waves do not rotate the polarization plane 
of the electromagnetic field, to first order ||. In this aspect, they behave like ordinary 
gravitational lenses. 
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4 Order of magnitude estimates 



In order to apply the previous theory, and for the multiple images to be detectable, the 
following conditions must be satified: 

1. the scale of separation 5 ~ h between different images must not be smaller than 
1CT 3 arcseconds, which gives 

/i>5-l(T 9 (4.1) 

at the deflection place; 

2. the impact parameter r must satisfy 

r > \ g . w . ; (4.2) 



3. the lens must not be exceptionally rare; 

4. the period of the lensing gravitational wave must not be too short (let us say 
P < 10 8 s), in order to appreciate variability of the images. 

Multiple imaging by gravitational waves is possible, and it occurs when the waves satisfy 
the approximate inequality (|1.7|) . We examine the validity of this condition for the most 
common astrophysical sources of gravitational radiation. 

Stellar core collapse: using the values of h and P predicted in studies of collapsing 
homogeneous ellipsoids, (|4.1| ) and ( fOj) allow only a rather narrow range for the impact 
parameter r, for which Dh/P ~ 10 if D ~ 6 ■ 10 15 cm. The late phase when the ellipsoid 
has settled down as a rapidly rotating neutron star does not give appreciable lensing. 
Asymmetry due to the core's bouncing gives Dh/P ~ 10 if D ~ 10 12 cm, r ~ 3 ■ 10 7 cm, 
or if D ~ 6 • 10 16 cm, r ~ 2 • 10 12 cm. Extrapolating the results obtained in studies of 
the perturbations of pressureless spherical collapse leading to the formation of a black 
hole, one obtains Dh/P ~ 10 if D ~ 3 ■ 10 16 cm (with a rather narrow range of values 
of r). 

Final decay of a neutron star/neutron star binary: rough estimates for the final decay 
of a binary system composed of two neutron stars give Dh/P ~ 10 if D ~ 3 • 10 16 cm, 
r ~ 5 • 10 10 cm, or if D ~ 6 • 10 7 cm, r ~ 6 • 10 7 cm. 

Black hole collisions: if two black holes collide with enough angular momentum to go 
into orbit before coalescing, one has Dh/P ~ 10 if D ~ 3 ■ 10 17 cm, r ~ 3 • 10 10 cm, or 
if D ~ 6 • 10 19 cm, r ~ 10 12 cm. 

The binary pulsar: the binary pulsar PSR 1913+16 is believed to radiate gravitational 



5 



waves in a continuous way with amplitude given, in order of magnitude, by h ~ Q/r ~ 
Ma 2 u 2 /r, where Q and M are the quadrupole moment and the mass of the pulsar, and a 
is the semimajor axis of the binary system. Conditions ( flip and ( |4.2[ ) are incompatible, 
hence multiple imaging is not possible in this case. 

The gravitational wave background: one has Dh/P « Jfl g . w .D/ R, where R is the radius 
of the universe and Q g . w . is the density of gravitational waves (in units of the critical 
density). Upper bounds on Q g . w . give fl g . w . < 1 for all frequencies, and Q g . w . <C 1 in 
many bands. Moreover, since D/R <C 1, also Dh/P <C 1, and multiple imaging by the 
gravitational wave background is, on average, impossible. 

As a conclusion, the creation of multiple images by a gravitational wave is possible 
and it is expected to occur. There are high amplification events associated to multiple 
images, in contrast with previous conclusions based on an inadequate formalism. The 
details of this work, and the study of the probability of observing multiple images and 
strong amplifications, and the details of the phenomenon will be the subject of a future 
publication |7| . 
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